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Abstract 

Using single cluster flip Monte Carlo simulations we accurately determine 

new finite size scaling functions which are expressed only in terms the variable 

X = £^l/L, where is the correlation length in a finite system of size L. Data 

for the d=2 and d=3 Ising models, taken at different temperatures and for 

different size lattices, show excellent data collapse over the entire range of 

scaling variable for susceptibility and correlation length. From these finite 

size scaling functions we can estimate critical temperatures and exponents 

with rather high accuracy even though data are not obtained extremely close 

to the critical point. The bulk values of the renormalized four-point coupling 

constant are accurately measured and show strong evidence for hyperscaling. 
PACS: 64.60, 05.70.Jk, 05.50.-Fq 
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I. INTRODUCTION 



True critical phenomena can possibly take place only in the limit in which the size of 
the system becomes infinite (i.e. the thermodynamic limit). The singular behavior of a sys- 
tem near a critical point is characterized by the bulk values of various physical quantities; 
it is technically impossible, however, to directly obtain information about infinite lattices 
from Monte Carlo simulation. Practically, however, it is not necessary to make the size 
of the lattice infinite in order to estimate a thermodynamic value through a Monte Carlo 
simulation of finite lattices: The concept of the finite-size-scaling (FSS) [|I|, introduced to 
extrapolate the information available from the finite system to the infinite volume limit, has 
been remarkably successful. The most frequent application of FSS has been primarily 
concerned with extracting some universal quantities such as the critical exponent v or some 
ratios of the critical exponents, without the knowledge of the bulk values in the scaling 
regime. The standard finite size scaling variable is x = th^^l^ with the reduced temperature 
t = \Kc — K\jKc^ where K is the inverse coupling K = J/kT and Kc is the inverse crit- 
ical coupling. Of course, use of this variable presupposes knowledge of the correct critical 
coupling and the uncertainty in Kc introduces a source of error into the analysis. Another 
formulation proposed by Fisher used t = {K^;' — K\/K^), but since this requires knowledge of 
multiple finite lattice "critical couplings" it has thus seldom been used. Some non-universal 
critical parameters like Kc can also be calculated by different FSS techniques, e.g. the 4*^ 
order cumulant ratio method P| or the microcanonical Monte Carlo method [Q. 

Nonetheless, determining bulk values (i.e. in the thermodynamic limit) is an important 
task of MC simulations, because physical quantities can then be directly compared with 
experiment. Also, the variation of a suitable thermodynamic quantity with temperature 
near criticality characterizes its critical behavior, even if it is not describable by a power 
law. A criterion telling whether a quantity measured on a finite lattice at a temperature T 
is distinguishable from the thermodynamic value (value in the thermodynamic limit) is the 
ratio of the linear size of the lattice (L) to the correlation length (,^(T)): provided L/^{T) 
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is sufficiently large, the measured quantity becomes essentially independent of L. Thus, 
one needs very large L at temperatures where ^(T) becomes large. Unfortunately, in this 
situation critical slowing down limits the quality of the data. Recently, new techniques of 
FSS have been introduced which enable one to extract correct thermodynamic values 
indirectly for a variety of physical quantities. The feature characteristic of these techniques 
is the calculation of some FSS functions defined in terms of a non-conventional FSS variable. 

In this paper, we numerically calculate certain FSS functions which are different from 
the "usual" ones and extract estimates for the values of critical parameters for the two 
and three dimensional Ising models. In the next section we provide theoretical background, 
and in the following section we calculate bulk values of the correlation length (^), magnetic 
susceptibility (x), and the renormalized four-point couphng constant {g^'^^)- We summarize 
and conclude in the final section. 

II. THEORY AND SIMULATION 

The fundamental assumption of FSS theory is that v4i(t), the value of some thermo- 
dynamic quantity A on a finite lattice of linear size L, can be expressed as 



for a bulk quantity A which has a power-law critical singularity A{t) ~ t^^ where t = 
{Kc — K)/Kc. Eq.(0) is valid for values of L and ^{t) which are large; otherwise, there 
should be corrections to FSS, which unless explicitly stated are ignored throughout this 
work. 

Notice that using the critical form for ^, ^{t) ~ t~^, we can rewrite the scaling variable 
s{L, t) as 



(1) 



s{L,t) = {A{t)/L^/n"'^'', 



(2) 



so that Eq.(0) may be rewritten as 
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AL{t) = A{t)J^A{s{L,t)), (3) 

where the relation between the scahng functions /a and Ta is given by 

:Fa{s) = s^/'^fAis). (4) 

For A = C,, Eq.(^) shows that C,L{t)/L is just a function of ^{t)/L and vice versa, and this 
leads to the relation 

AL{t) = A{t)QA{xmt)), (5) 

where x{L,t) = ^L{t)/L is the ratio of the correlation length on a finite lattice to the 
lattice size, and Qa{,x) is given by 

QA{x)=TAUi\^))- (6) 

Using the same observation, it is trivial to obtain |^ another equivalent form, 

A^L{t)=AL{t)GA{x{L,t)), (7) 

where 6 is a scaling factor and Qa{x) is another scaling function. 

It is evident that given one can determine the other two scaling functions from Eqs.(p 
and d^), and all the scaling functions, /a, J^a, Qa, and Qa should be universal. It has also 
been argued that a certain asymptotic form of fA{s) can be expressed in terms of the 
critical exponent 5; by fitting this functional form one can extract an estimate for the critical 
exponent. 

It is worth stressing that use of the scaling function Q rather than would be more 
convenient in many cases, particularly because one does not need the bulk correlation length 
to define the former. Note that there is no explicit t dependence of the scaling variables, 
so that knowledge of the critical coupling is not required, and that x becomes independent 
of L at criticality. This L independent value of x at criticality, Xc-, which characterizes a 
universality class for a given geometry [^, forms the upper bound of x. In other words, 
the scaling function Q is defined only over < x < Xc. A priori, the two limits of the 
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scaling function Q are known for a continuous phase transition: lima;^o Q{x) 



1 and 



lima-^a;^ Q ^ 0, because Al converges to its bulk value in the former case while A{t) diverges 
in the latter case with AL{t) finite. In general, as we will show in this work, for A = ^, x, 
or g^^\ Qa{x) turns out to be a monotonically decreasing function of x. 

It is important to realize that the knowledge of the scaling function Q near x ~ plays 
as relevant role as that near x ~ Xc to the extraction of necessary information of the critical 
behavior in (deep) scaling region. It can be easily seen by noting that x{L, t) for a fixed 
temperature arbitrarily close to criticality can be made arbitrarily close to zero by simply 
choosing a value of L sufficiently large. 

Eqs.(|D and (j^) do not include any critical exponents, so that one might conjecture that 
their validity can be extended to non-power-law singularities. Although a general proof 
of this conjecture is missing, Liischer obtained an explicit expression for the inverse 
correlation length (mass gap), that is consistent with Eq.(^), for the two-dimensional (2D) 
0{N) {N > 2) spin models which exhibits an exponential critical singularity. Also, very 
extensive numerical verification |^ of Eq.(^ was given for x and ^ for the 2D 0{N) models 
with N=2 and 3. 

In order to define a correlation length, we consider the Fourier transform of the (con- 
nected) two-point correlation function. 



where S'x denotes the spin variable at site x. When x is sufficiently large, < • S'x >c~ 



G(k) = ^ exp{ik ■ x) < Sq ■ S:^. > 



(8) 



X 



e l^l/^i- holds [[T^], so we will have 



G(k) 



Gior'[i+k'eL+oik')] 



(9) 



By choosing k = (27r/L, 0), we obtain 




(10) 



for values of L that are large enough that terms ~ fc^ 



can be safely ignored. 
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The renormalized coupling constant gf^ may be defined as [12 



gf = 3(L/a)^ f/L, (11) 

where D is the lattice dimensionality and the 4th order cumulant is given by f/^ = 1— < 
> /3 < 5*^ with S being the order parameter. The bulk g^^^ has a well defined scaling 
behavior 

g^^^ describes the non-Gaussian character of the the model, i.e., only for a Gaussian model 
does g^^\t) vanish as t — > in the absence of certain multiplicative correction to scaling, 
implying the violation of the hyperscaling relation Du — 2A + 7 > 0. For a system where 
the hyperscaling relation is satisfied (without certain multiplicative correction to scaling as 
in the four dimensional Ising model), its bulk value in the scaling regime remains a constant 
that characterizes its universality class. 

Employing the single cluster Monte Carlo algorithm we simulated the 2D and 3D 
Ising models, on the square and simple cubic lattices respectively, with fully periodic bound- 
ary condition imposed. For each lattice at a given temperature, we generated up to 30 bins 
of data each of which is composed of 10000 measurements. In order to reduce the correla- 
tion between data points, only configurations 3-7 Monte Carlo step apart were considered. 
Our quoted errors, which are purely statistical in nature, are the standard deviation of the 
binned values. Aware of the bad performance of some random number generators in the 
context of the single cluster algorithm |1^, we have double-checked our results by com- 



paring data generated by two different implementations of the algorithm each one using 
a different kind of random number generator. Most of data were obtained with a linear 
congruential random number generator of the form Xj+i = 69069 x + 1 mod (2^^). The 
other random number generator was a multiplicative, lagged Fibonacci generator of the form 
Xj+i = Xj_4423 X Xj_i393, which showcd a good performance in a single cluster simulational 



test of Ising system [0. We observed complete agreement between the two sets of data 



within our statistical errors. We therefore, beheve that to within error bars quoted here our 
data are not biased due to any correlation among the random numbers. (We also tried the 
well known R250 routine but the data exhibited some systematic deviation and we did not 
consider them in our analysis.) 



III. RESULTS 
A. 2D Ising Model 

We now investigate the finite size behavior for a variety of multiplicatively renormal- 



izable physical quantities |]I0| defined on a finite lattice of linear size L, in particular the 
susceptibility x ^i-nd correlation length ^. First, we choose a certain K and perform mea- 
surements of Al{K) for various values of L. In Table(l) we present our data for K = 0.425, 
with L varying from 16 to 150. The reason for starting from L = 16 is that near this L the 
systematic error in our definition of Eq-(0) is about ~ 10~^ which is comparable to our 
typical statistical errors. From these data in Table(l) one sees that ^ is indistinguishable 
from bulk value for L > 80, to within very small statistical error. In terms of the scaling 
variables, this means that s > 5.076 or a; < 0.196. As we stressed earlier, this condition 
holds (in terms of the scaling variables) regardless of the temperature as long as it remains 
in the scaling regime. (This is indeed the fundamental statement of FSS.) We note that for 
X the thermodynamic condition holds for a slightly larger value of s (smaller x) than that 
for ^, namely s > 6.346. Figure(l) shows the data-collapse for J-'a{s). 

From the data in Table(l), we can easily determine J-'a{s) and Qa{x). In order to satisfy 
the asymptotic conditions, for the former one may try either a polynomial function of 1/s 
or e~^, while for the latter a polynomial of x or e~^^^ may be tried. That is, 

J^a{s) = 1 + bi/s + + . . . , (13) 
Qa{x) = l + bix + hx^ + . . . , (14) 

or 



J^a{s) ^ 1 + cie-' + C2e-^' + . . . , (15) 
Qa{x) = 1 + cie-^/" + C2e-2/" + . . . . (16) 

In general, it turns out that for the same number of fitting parameters the polynomial 
of the exponentials fits better than that of the simple scaling variables. This is espe- 
cially true for the magnetic susceptibility and the four-point renormalized coupling con- 
stant. For instance, by considering terms up to the fourth order of the polynomial, 
X^/-^DF(degree of freedom) = 3.3 and 0.3 respectively for the Qxi^) assuming the sim- 
ple polynomial and that of the exponential. Considering up to the e~^/^ term, we obtain 
ci = -2.402, C2 = -16.338, C3 = 80.688, C4 = -134.6 for with x^/Ndf^0.31, while they 
are -0.768, -8.490,31.032, and -89.203 respectively for with x^/Ndf=0-20. 

A priori, the estimates are accurate only for x < ^Lq/Lq = xq, with Lq denoting the 
smallest value of L ior K — 0.425; for the estimate of the coefficients for x > Xq, one 
needs similar data for Al at a. larger K, which might modify the values of the coefficients. 
Nevertheless, with the information of the finite size-scaling function Qa, it is now possible to 
extract accurate bulk data from the Monte Carlo data on a modestly small lattice provided a 
data point (at another temperature) satisfies x < xq. Our results for K = 0.430, 0.434, 0.436, 
and 0.438 are summarized in Table(2). 

We notice that the bulk values thus extracted for a given K do not change with respect 
to L, indeed verifying this form of FSS for the model (see Figure(2) also). We also note that 
the values of the bulk ^ thus extracted are in excellent agreement with the corresponding 
exact values given by the formula 

1/^(K) = In(cothX) - 2K, (17) 

within typically less than 0.5% of the statistical errors. As a test of our x, we fitted the data 
over the range from K — 0.425 to 0.438 to 

X ~ t-\ (18) 
The best fit gives Kc = 0.44070(5) and 7 = 1.755(9) (7=1.752 by fixing Kc to the exact 
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critical point), being extremely close to the exact values. Since all of our data used for the 
analysis were for t > 0.006, the quantity of the result is surprisingly good. 



B. The 3D Ising Model 

We begin with our Monte Carlo measurement at K = 0.220, the results of which are 
summarized in Table(3). We observe that the thermodynamic condition for the C, is almost 
satisfied for L/^l ^ 60/10.89 ~ 5.5. Nevertheless, we also note that x ^^(^ 9^'^^ increase, 
albeit very slowly, beyond this value. This is another indication that certain quantities 
converge to the thermodynamic value more slowly than the correlation length. 

Assuming that Al{K = 0.220) reaches its bulk value for L = 70, we obtain 

ci ~ -0.418, C2 ~ -18.83, C3 ~ 99.38, C4 ~ -436.4, (19) 
ci ~ -0.607, C2 ~ -56.31, C3 ~ 416.75, C4 ~ -1399.9 (20) 

ci ~ -7.238, C2 ~ 21.42, C3 ~ -16.09, C4 ~ 3.67, (21) 

respectively for Q^{x), Q^{x), and Qg(4)(x), over the range < x < xq = 0.491. 

Based on the knowledge of Qa{x), we calculate the bulk values of the correlation length, 
magnetic susceptibility, and four-point renormalized coupling constant for various K up to 
K = 0.2212. The largest value of L we simulated for the calculation is just 64. Obviously, 
the computation of the bulk value at a larger K (than 0.2212) requires a larger value of L 
in order to keep the value of x smaller than xq- One way to avoid the need for a larger 
L is to repeat the measurement of at a slightly larger value of K, e.g. K = 0.2212; 
this will extend the range of x over which Qa is accurately computed. In the region where 
Qa{x) ~ 0, however, one needs extremely precise measurements to reduce the errors in the 



estimates of the bulk values p6[|. 



A summary of our results is shown in Table(4). We see that the four- point renormalized 
coupling constant remains unchanged, i.e., g^^^ ~ 24.5, for K > 0.2206. Its slow variation 
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for K < 0.2206 may be due to the presence of correction to scaling. Hence the hyperscahng 
relation is indeed satisfied, confirming the previous verification (within rather larger statis- 
tical errors though) based on the traditional Monte-Carlo measurement ||12| We would 
like to stress, however, that it was not possible to measure g^^^ beyond K = 0.2206 in Ref. 



T^, even using L as large as L = 90. The bulk x and ^ thus extracted are compared with 
those traditionally obtained, again yielding remarkable agreement (see Table(4)). Figure(3) 
exhibits excellent data-collapse for the finite size scaling function Qa{x) for A = x, ^, and 

In order to determine Kc, ly, and 7 we fitted our bulk data over the range 0.217 < 
K < 0.2212 to the simple power-law singularity. We fixed the critical point in the fit, 
and then repeated the fit for several different fixed critical points. The results, shown 
in Fig. 4, indicate that the values of the ^ and x data favor the range of over 
0.221640 < Kc < 0.221670, being consistent with recent other results |1^. The empirical 
formulae we obtained from the best fit are as follows: 

^ = C^{\K -Kc\/Kcy, Q = 0.4710, 7^^ = 0.221658, z/ = 0.6418 (22) 
X = C^i\K - Kc\/Kc)~^, = 1.0892, = 0.221646, 7 = 1.2388 (23) 

The value of u is larger by approximately 1.5 percent than those extracted by most other 
methods, while the value of 7 agrees up to ~ 10^^. The effect of including the term of the 
confluent correction to the scaling turns out to be minimal: the confluent correction-term 
would usually be important when data with rather smaller bulk values of the correlation 
length are considered. Given the modest size lattices used and the distance from the critical 
points at which the actual measurements were made (t > 0.002), we flnd that the agreement 
with high resolution studies to be extremely gratifying. 



IV. DISCUSSION AND CONCLUSION 



In this paper we computed an alternative flnite size scaling function, deflned in terms of 
the scaling variable C,l/L, for the 2D and 3D Ising models. This type of flnite size scaling 
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function has the advantage of being defined even for small values of L, without any prior 
information on the critical behavior of the system. Thus, our procedure can test FSS itself 
by means of data collapse as far as can be accurately measured; in this manner, we 
observed that the effect of the violation of FSS is negligibly small at least for L > 16 for the 
two and three dimensional Ising models. We illustrated how the function can be used for 
the extraction of correct bulk values near criticality, and that it can be used in extracting 
accurate critical parameters provided the values of the correlation length are sufficiently 
large, i.e. approximately ^ > 5. 

One might wonder if this technique requires an accurate bulk value of a physical quantity 
A at least at a point of temperature. As far as the extraction of critical parameters is 
concerned, however, this is not necessarily the case. To see this, imagine that we start with a 
fake 'bulk value' A'{t) instead of correct one A{t). The scaling function Q' defined in terms of 
the fake bulk value, Q' a = AL{t) /A'{t) = -^^Qa, simply rescales the correct scaling function 
by a constant; accordingly, every bulk value calculated at any other temperature using Q'^ 
rescales the correct one with an overall factor '^^^^ overall factor is unimportant for 

the extraction of the critical behavior. One can thus repeat our analysis arbitrarily close to a 
critical point, where the effect of correction to scaling can be arbitrarily small. We anticipate 
that such an analysis will yield extremely accurate estimate of the critical parameters, and 
our study on the 3D Ising model along this line is under way. (This observation is also 
important for some calculations of lattice gauge theory. For example, in full lattice QCD, 
the computation of the ratio of the mass of various lattice 'hadrons' is of primary concern, 
and for this purpose the overall factor is simply unimportant.) 

We would like to stress again that the technique we have illustrated is extremely general; 
it holds regardless of the functional form of the critical singularity, and irrespective of the 
quantity as far as it is multiplicatively renormalizable. We demonstrated this point taking 
the example of the four-point renormalized coupling constant, whose bulk value is notori- 
ously difficult to measure through traditional Monte Carlo simulation |]T^ ||T3. Although 



our current estimate of the critical parameters cannot compete with the highest resolution 
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Monte Carlo studies, our estimates are really surprisingly good considering how far from 
the critical point the data are taken. Thus, although the determination of a finite lattice 
correlation function is needed, this method offers a simple alternative to standard finite size 
scahng methods. 
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TABLES 

TABLE I. Al as a function of L at ^-=0.425. ^oo(A'=0.425)=15.7582... from Eq.(17). Note 
that our converges to its bulk value for L > 80, within the small statistical errors, and that the 
value of X monotonically decreases with L. 



L 




X 


XL 


16. 


9.83(3) 


0.614(2) 


102.4(2) 


18. 


10.56(4) 


0.587(2) 


119.8(3) 


20. 


11.18(5) 


0.559(2) 


137.0(6) 


22. 


11.74(4) 


0.534(2) 


153.5(5) 


25. 


12.44(5) 


0.498(2) 


177.1(7) 


27. 


12.85(5) 


0.476(2) 


191.9(7) 


30. 


13.37(6) 


0.446(2) 


212.7(6) 


32. 


13.69(7) 


0.428(2) 


224.7(9) 


34. 


13.92(6) 


0.409(2) 


235.2(9) 


36. 


14.19(6) 


0.394(2) 


246.6(8) 


40. 


14.54(6) 


0.363(2) 


264.9(6) 


50. 


15.19(8) 


0.304(2) 


296.0(1.0) 


60. 


15.40(6) 


0.257(1) 


312.5(1.0) 


70. 


15.62(7) 


0.223(1) 


321.3(1.4) 


80. 


15.71(8) 


0.196(1) 


326.6(1.3) 


100. 


15.75(10) 


0.157(1) 


329.4(1.8) 


110. 


15.77(10) 


0.143(1) 


331.0(1.5) 


150. 


15.75(14) 


0.105(1) 


331.4(1.6) 
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TABLE II. The extraction of the ^ and x for the 2D Ising model based on the computation 
of the Qa{x)- The 'ave.' for each K denotes the average over the extracted bulk values from the 
different values of L{x). 



K 


L 




X 


XL 




X 


0.430 


20. 


13.00(7) 


.650 


158.5(4) 


23.16(17) 


650.4(16.2) 




30. 


16.67(7) 


.556 


274.5(4) 


23.33(11) 


653.2(5.4) 




40. 


18.97(8) 


.474 


377.3(6) 


23.20(9) 


647.5(4.2) 




50. 


20.57(8) 


All 


461.1(8) 


23.31(8) 


650.2(3.8) 




oO. 


21.57(8) 


.359 


524.5(1.5) 


o o on / ''7\ 

23.30(7) 


651.0(5.5) 




80. 


22.57(10) 


.282 


597.0(2.2) 


23.23(8) 


649.7(6.3) 




100. 


22.89(11) 


.229 


625.3(2.3) 


23.15(9) 


646.6(5.9) 




120. 


23.12(17) 


.193 


638.2(5.0) 


23.23(13) 


647.2(12.4) 


ave. 


- 


- 


- 


- 


23.21(11) 


649.5(2.2) 


exact 










23. 22 




0.434 


80. 


32.94(7) 


.412 


1047.5(2.5) 


37.34(7) 


1478.6(11.9) 




160. 


36.77(9) 


.230 


1426.5(9.0) 


37.19(14) 


1476.0(23.1) 


ave. 


- 


- 


- 


- 


37.27(11) 


1477.3(1.8) 


exact 


- 


- 


- 


- 


37.21 


- 


0.436 


50. 


31.58(9) 


.632 


761.7(1.4) 


53.10(20) 


2745.0(43.7) 




60. 


35.35(11) 


.589 


977.4(2.1) 


53.18(19) 


2743.8(39.8) 




70. 


38.54(12) 


.551 


1186(3) 


53.41(18) 


2759.0(37.7) 




80. 


41.23(13) 


.515 


1389(3) 


53.65(17) 


2779.4(31.8) 




90. 


43.08(7) 


.479 


1563(2) 


53.01(8) 


2724.8(17.5) 




100. 


45.04(13) 


.450 


1746(2) 


53.37(14) 


2751.1(24.7) 




120. 


47.96(14) 


.400 


2033(5) 


53.69(13) 


2776.8(22.4) 
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160. 


50.55(20) 


.316 


2384(7) 


52.96(17) 


2724.1(22.9) 


ave. 










53.30(28) 


2750.6(20.8) 


exact 










53.16 




0.438 


80. 


51.88(13) 


.649 


1777(4) 


91.97(31) 


7209.7(142.4) 




160. 


76.31(34) 


.477 


4259(13) 


93.66(39) 


7378.2(93.0) 


ave. 










92.82(1.20) 


7294.0(119.1) 


exact 










92.86 




TABLE III. AiiK) at K = 0.220 for the 3D Ising model. 


L 








X 


XL 


(4) 

9l 


16. 




7.85(2) 




0.491(1) 


228.8(7) 


9.52(5) 


20. 




8.85(2) 




0.443(1) 


298.1(7) 


11.43(5) 


24. 




9.56(2) 




0.398(1) 


351.8(1.1) 


13.44(7) 


30. 




10.20(3) 




0.340(1) 


407.1(1.1) 


16.5(1) 


36. 




10.56(2) 




0.293(1) 


439.5(1.2) 


19.2(2) 


40. 




10.68(3) 




0.267(1) 


455.2(1.5) 


21.4(3) 


50. 




10.85(3) 




0.217(1) 


467.9(1.3) 


23.5(5) 


60. 




10.89(3) 




0.182(1) 


472.3(0.8) 


24.2(5) 


70. 




10.91(3) 




0.156(1) 


473.0(1.1) 


24.7(5) 
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TABLE IV. The extraction of the ^ and % for the 3D Ising model. The uncertainties in the 



values of x are not considered for our error estimates. Traditional Monte Carlo measurements are 
in rows labeled MC. 



K 


L 




X 


XL 


(4) 

9l 




X 


9^'^ 


0.217 


16. 


5.30(1) 


.331 


114.6(2) 


17.7(1) 


5.62(1) 


130.8(2) 


25.5(1) 




20. 


5.49(1) 


.275 


124.5(4) 


21.2(2) 


5.63(1) 


130.7(4) 


25.4(1) 




28. 


5.60(1) 


.200 


130.1(1) 


24.6(2) 


5.62(1) 


131.0(3) 


25.8(2) 


ave. 












5.62(1) 


130.9(4) 


25.6(3) 


0.219 


20. 


7.33(2) 


.366 


212.4(7) 


15.3(2) 


8.04(2) 


262.2(9) 


24.9(6) 




30. 


7.88(3) 


.263 


251.5(1.5) 


21.8(5) 


8.03(4) 


260.9(1.6) 


25.7(6) 


ave. 












8.03(4) 


261.6(2.1) 


25.3(1.1) 


MC 


50 


- 


- 


- 


- 


8.1 


263 


26 


0.2203 


30. 


11.16(8) 


.372 


481.4(5.3) 


14.8(1) 


12.31(9) 


602.7(6.6) 


24.6(2) 




40. 


12.00(3) 


.300 


561.4(2.6) 


19.3(2) 


12.44(3) 


608.2(3.1) 


25.1(4) 


ave. 












12.38(15) 


605.5(8.8) 


24.9(7) 


0.2206 


24. 


11.20(5) 


.467 


467.4(2.8) 


10.3(2) 


14.55(6) 


834.2(5.0) 


24.5(2) 




30. 


12.49(10) 


.416 


589.2(2.7) 


12.5(2) 


14.66(12) 


844.4(3.9) 


24.6(2) 




36. 


13.28(5) 


.369 


673.8(3.6) 


14.9(3) 


14.60(5) 


836.8(4.5) 


24.5(3) 




40. 


13.70(4) 


.343 


718.0(3.4) 


16.3(2) 


14.66(4) 


838.4(5.0) 


24.4(3) 




50. 


14.28(5) 


.286 


789.7(4.0) 


19.4(3) 


14.68(5) 


839.7(5.3) 


24.2(4) 




60. 


14.38(3) 


.240 


809.9(1.6) 


22.4(3) 


14.53(5) 


827.6(2.6) 


25.1(4) 


ave. 












14.61(6) 


837.0(5.7) 


24.6(3) 


MC 


75 










14.5 


828 


24 
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0.2210 40. 16.75(3) 

50. 18.15(4) 

ave. 

0.2212 56. 21.98(7) 

64. 22.89(7) 

ave. 



.419 1045.0(3.1) 

.363 1245.7(4.1) 

.393 1796.0(10.5) 

.358 1971.6(10.4) 



12.3(1) 19.73(5) 
15.1(2) 19.82(8) 
19.77(12) 

13.6(1) 24.89(10) 
15.5(2) 24.86(10) 
24.88(12) 



1511.5(5.5) 24.5(2) 

1524.3(7.1) 24.4(4) 

1518.(15.) 24.5(3) 

2382(14) 24.5(2) 

2382(13) 24.4(4) 

2382(14) 24.5(3) 



Figure Captions 

Fig.(l): !Fa{s) for the 2D Ising model. Each 'curve' demonstrates the data collapse for 
two different values of K. Note that the lower curve converges to Ta{s)—'^, i.e., thermody- 
namic hmit more slowly than the upper curve. 
Fig. (2): Qa{x) for the 2D Ising model. 

Fig. (3): Qa{x) for the 3D Ising model. The solid hues are from the best fits. 
Fig. (4): The resulting estimates for the critical exponents 7 and v for different choices 
of the critical coupling Kc- The correspondent errors are given by the light dotted lines. 
Solid (dashed) curves correspond to x (^) data. The shadings show acceptable values for 
the critical parameters. On the bottom: plotted against Kc] vertical arrows locate our 
best estimates for K^i whereas the horizontal ones indicate the error bars. 
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